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Abstract 

The ground-state phase diagram of the spin-one- half Falicov-Kimball model (FKM) 
is studied in the one dimension using small-cluster exact-diagonalization calculations. 
The resultant exact solutions are used to examine possibilities for valence and metal- 
insulator transitions in this model. A number of remarkable results are found, (i) 
The phase separation in the spin-one-half FKM takes place for a wide range of /- 
electron concentrations rif and d-f interactions G, including G large, (ii) In the strong 
coupling limit (G > 4) the model exhibits a pressure induced discontinuous insulator- 
metal transition from an integer- valence state (n/ = 1) into another integer- valence 
state (rif = 0). (iii) For intermediate values of G (G ~ 2.5) the FKM undergoes a 
few discontinuous intermediate-valence transitions. There are several discontinuous 
insulator-insulator transitions from rif = 1 tony = 1/2 and a discontinuous insulator- 
metal transition from rif = 1/2 to rif =0. (iv) In the weak coupling limit (G < 2) the 
model undergoes a few consecutive discontinuous and continuous intermediate- valence 
transitions as well as a discontinuous metal-insulator transition. 

PACS nrs.:75.10.Lp, 71.27.+a, 71.28.+d, 71.30.+h 



1 



1 Introduction 



Recent new measurements of the insulator- metal transitions in SmBg jlj] and in the transi- 
tion-metal halides [|J, have once again sparked an interest in valence and insulator- metal 
transitions |J. These transitions are observed in a wide group of substances formed by 
transition-metal oxides as well as rare-earth sulfides and borides, when some external pa- 
rameters (like pressure or temperature) are varied. They are in many cases first-order 
phase transitions, however, second-order transitions ranging from very gradual to rather 
steep are also observed 0]. 

To describe all such transitions in a unified picture Falicov and Kimball introduced 
a simple model in which only two relevant single-electron states are taken into account: 
extended Bloch waves and a set of localized states centered at the sites of the metallic 
ions in the crystal. It is assumed that insulator-metal transitions result from a change 
in the occupation numbers of these electronic states, which remain themselves basically 
unchanged in their character. The Hamiltonian of the model can be written as the sum of 
four terms: 



ijtr iaa' 

Ef fiaficr + ^" E f^rficrfiLufi-a, (1) 
icr 2(7 

where fi a are the creation and annihilation operators for an electron of spin a in the 
localized state at lattice site i with binding energy Ef and df a , d{ a are the creation and 
annihilation operators of the itinerant electrons in the d-band Wannier state at site i. 

The first term of (1) is the kinetic energy corresponding to quantum-mechanical hopping 
of the itinerant d electrons between sites i and j. These intersite hopping transitions are 
described by the matrix elements tij, which are —t if i and j are the nearest neighbors 
and zero otherwise (in the following all parameters are measured in units of t). The 
second term represents the on-site Coulomb interaction between the d-band electrons with 
density = N^/L = J^ia df a di a and the localized / electrons with density nf = Nf/L = 
t J2ia fiafio, where L is the number of lattice sites. The third term stands for the localized 
/ electrons whose sharp energy level is Ef. The last term represents the intra-atomic 
Coulomb interaction between the localized / electrons. 
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In spite of the fact that the FKM is one of the simplest examples of interacting fermionic 
system, the theoretical picture of valence and insulator-metal transitions remains still un- 
certain in the framework of this model. Even, in the existing literature on this model, 
different answers can be found on the fundamental question whether the FKM can de- 
scribe both the discontinuous and continuous changes of the / (d)-electron occupation 
number nj (n^) as a function of the /-level energy Ef §. It should be noted that this 
question is indeed crucial for the systems mentioned above, since, supposing [|J] that the ex- 
ternal pressure shifts the energy level Ef, the valence changes observed in some rare-earth 
and transition-metal compounds (SmS, SmBg, Ti203, and so on) could be understandable 
purely electronically. Unfortunately, it was found that valence and insulator-metal transi- 
tions are very sensitive to the approximation used. Various approximations || (mean-field, 
virtual crystal, CPA, etc.) yield very different and often fully controversial results. This 
indicates that the study of valence and insulator-metal transitions may be successful only 
with methods which are relatively insensitive to the type of approximation used and, of 
course, with exact methods. 

In our previous papers |7|, || we have showed that the method of small-cluster exact- 
diagonalization calculations is very effective in describing ground-state properties of both 
the spinless and spin-one-half FKM. For the spinless FKM H the exact numerical calcu- 
lations (over the full set of /-electron configurations) can be performed on relatively large 
clusters (L ~ 36), even without some special computational effort. For such clusters the 
finite-size effects are practically negligible and the results can be satisfactory extrapolated 
to the thermodynamics limit (L — > oo). Using this method we have described success- 
fully the ground-state phase diagram as well as the picture of valence and metal-insulator 
transitions in the spinless FKM |jj . The situation for the spin-one-half FKM is more com- 
plicated. The full Hilbert space of the spin-one-half FKM is much larger than one of the 
spinless FKM and this fact impose severe restrictions on the size of clusters that can be 
studied by the exact-diagonalization method. Our recent numerical computations per- 
formed for the spin-one-half FKM with U and G finite showed that clusters with L > 24 
are beyond the reach of present day computers. Fortunately, the size of the Hilbert space 
can be reduced considerably in some special, but physically still interesting limits, e.g., 
U — > oo. In the limit U — ► oo states with two / electrons at the same site are projected out 
thereby much larger clusters (L ~ 36) become accessible for the numerical investigation in 
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this reduced subspace. For this reason all calculations presented in this paper have been 
done at U = 00. The main goal for performing these calculations was to construct the 
comprehensive phase diagram of the spin-one-half FKM. The second goal of our numerical 
study was to find and describe all possible types of valence and metal-insulator transitions 
in this model. We show that both the ground-state phase diagram as well as the picture 
of valence and metal-insulator transitions obtained for U = 00 strongly differ from ones 
obtained in our previous paper || for a restricted set of /-electron configurations and U 
finite. In particular, the exhaustive numerical studies of the model (at U = 00) performed 
on finite clusters up to 36 sites revealed some new unexpected features like the phase 
separation for all nonzero G and discontinuous metal-insulator transitions for G small, to 
mention only a few. In the present paper we discuss these features in detail. 

Since the /-electron density operators ff a fm of each site i commute with the Hamilto- 
nian (1), the /-electron occupation number is a good quantum number, taking only two 
values, Wi = 0, 1 according to whether the site i is unoccupied or occupied by the localized 
/ electron (configurations with W{ = 2 are projected out due to U = 00). Therefore the 
Hamiltonian (1) can be replaced by 

H = E hijdtdje + E f ^2wi, (2) 

ija i 

where h{j = Uj + GwiSij. 

Thus for a given /-electron configuration w = {wi,W2 ■ ■ ■ Wl} defined on a one-di- 
mensional lattice with periodic boundary conditions, the Hamiltonian (2) is the second- 
quantized version of the single-particle Hamiltonian h(w) = T+GW, so the investigation of 
the model (2) is reduced to the investigation of the spectrum of h for different configurations 
of / electrons. Since the d electrons do not interact among themselves, the numerical 
calculations precede directly in the following steps (we consider only the case Nf + = 
L, which is the point of the special interest for valence and metal-insulator transitions 
caused by promotion of electrons from localized / orbitals (f n — > / n_1 ) to the conduction 
band states), (i) Having G, Ef and w = {wi,W2 ■ ■ ■ Wl} fixed, find all eigenvalues 
of h(w) = T + GW. (ii) For a given Nf = J2i w i determine the ground-state energy 
E{w) = Ylk=i S + EfNf of a particular /-electron configuration w by filling in the 
lowest Nd = L — Nf one-electron levels (the spin degeneracy must be taken into account), 
(iii) Find the w° for which E(w) has a minimum. Repeating this procedure for different 
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values of Nf, G and Et , one can immediately study the ground-state phase-diagram of the 
model in the Nf-G plane as well as the dependence of the /-electron occupation number 
Nj = J2i w i on * ne /-level position Ef (valence transitions). 

2 Phase separation 

Let us first discuss the ground-state phase-diagram of the spin-one-half FKM. To reveal 
the basic structure of the phase diagram in the Nf-G plane {Ef = 0) we have performed 
an exhaustive study of the model on finite (even) clusters up to 36 sites. For fixed L 
the numerical calculations have been done along the lines discussed above with a step 
ANf = 2 and AG = 0.05. The results of numerical computations are summarized in 
Fig. 1. These results show that the phase diagram of the spin-one-half FKM consists of 
three main domains: the most homogeneous domain MHD (in Fig. 1 denoted as ■) and two 
phase separation domains PSD\ (denoted as o, x and +) and PSD2 (denoted as A). In 
the MHD the ground states are configurations in which the atomic or n-molecule clusters 
of /-electrons are distributed in such a manner that the distances between two consecutive 
clusters are either d or d + 2. Furthermore the distribution of the distances of d and d + 2 
has to be most homogeneous. Two basic types of the ground state configurations that 
fill up practically the whole MHD are displayed in Table. 1. In the PSD\ the ground 
state configurations are mixtures w&Lw e of aperiodic (periodic) configurations w and the 
empty configuration w e = {00 . . . 0} (see Ref. [|| for a definition of mixtures), i.e., all 
/-electrons are distributed only in one part of the lattice (configuration w) while another 
(connecting) part is free of /-electrons (phase separation). We have found three basic types 
of configurations w which form these mixtures: (i) the aperiodic atomic configurations w a = 
{lOfcjlOfcj . . . lOj^l . . . OfcjlOfol} (with ki > 0), (ii) the aperiodic nj-molecule configurations 
Wb = {lni0fc 1 l n2 0fc 2 • • • l ni fej l n . . . .0 fc2 l„ 2 fel l ni } (with 1 < rii < N f /2 and h > 0), and 
(hi) the iVj-molecule (segregated) configurations w c = {11 ... 1}. Three different regions 
of stability corresponding to mixtures w a &iw e , Wb&zw e and w c &iw e are denoted in Fig. 1 
as o, x and +. It is seen that the mixtures of the atomic configuration w a and the empty 
configuration are stable only at low /-electron concentrations and the Coulomb interactions 
G < 2.2. A direct comparison of results obtained for the spin-one-half and spinless FKM || 
shows that this region corresponds roughly to a region of phase separation in the spinless 
FKM. Outside these regions the phase diagrams of the spin-one-half and spinless FKM are, 
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however, strongly different. While the phase separation in the spinless FKM takes place 
only for weak interactions (G < 1.2) the spin-one-half FKM exhibits the phase separation 
for all Coulomb interactions. Even with increasing G the phase separation shifts to higher 
/-electron concentrations. Particularly, in the region 2.5 < G < 2.7 where the ground 
states are the mixtures of the n-molecule configurations Wb with the empty configuration 
the phase separation takes place for all Nf < L/2 and in the region G > 2.7 where the 
ground states are the segregated configurations ws = w c $zw e even for all Nf < L. At 
large /-electron concentrations n/ but in the opposite limit (G < 0.4) there exists another 
small domain of the phase separation PSD2 (denoted as A). The numerical results on 
finite lattices up to 36 sites revealed only one type of configurations that can be the ground 
state configurations in this domain, and namely, the mixtures of the periodic n-molecule 
configurations Wd with the fully occupied lattice u>/ = {ll...l} (the length of a connected 
cluster of occupied sites in these mixtures is at least L/2). 

The second step in our numerical studies has been the extrapolation of small-cluster 
exact diagonalization results on large lattices. In Fig. 2 we present the ground state phase 
diagram of the spin-one-half FKM obtained for L = 100 on the extrapolated set of con- 
figurations that includes practically all possible types of the ground-state configurations 
found on finite lattices up to 36 sites. In particular, we have considered (i) the most homo- 
geneous configurations Wh of the type a and b (see Table 1), (ii) all mixtures w a hw e with 
ki smaller than 6, (iii) all mixtures Wb&zw e with and k{ smaller than 6, (iv) all mixtures 
WdSzwf with periods smaller than 12, and (v) all segregated configurations. One can see 
that all fundamental features of the phase diagram found on small lattices hold on much 
larger lattices too. Of course, the phase boundaries of different regions corresponding to 
mixtures w a &zw e (a), Wb&w e (b), w c k,w e (c) and Wd&Wf (d) are now more obvious. Since 
the mixtures w a Szw e , WbSzw e , w c &iw e and Wd&Wf are metallic Q one can expect that the 
phase boundary between the MHD and PSD is also the boundary of the correlation induced 
metal-insulator transition. To confirm this conjecture it is necessary to show that the most 
homogeneous configurations from the MHD are insulating, or by other words that there is 
a finite energy gap A at the Fermi energy in the spectra of these configurations |1C]. The 
numerical results for the most-homogeneous configurations of the type a and b that fill 
up practically the whole MHD are displayed in Fig. 3. The insulating character of these 
configurations is clearly demonstrated by the finite A that exists for all nonzero /-electron 
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concentrations rif and Coulomb interactions G. Thus we can conclude that the spin-one- 
half FKM undergoes (on the phase boundary between the MHD and PSD) the correlation 
induced metal-insulator transition that is accompanied by a discontinuous change of the 
energy gap A. 

3 Valence and metal-insulator transitions 

The existence of a large metallic domain in the strong coupling region is the main difference 
between the phase diagram of the spin-one-half and spinless FKM. In the spinless FKM M 
the existence of the metallic phase was restricted only on a small region G < 1.2 and 
rif < 1/4 (rif > 3/4), while the remaining part of the phase diagram was insulating. The 
phase diagram of the spin-one-half FKM has a more complicated structure. In addition to 
the insulating phase corresponding to the homogeneous configurations Wh there are four 
metallic phases corresponding to four different classes of the ground state configurations: 
w a Szw e , Wb&w e , w c hw e and Wd&zwf. Of course, this fact has to lead also to different picture 
of valence and metal-insulator transitions induced by pressure (increasing Ef). Let us now 
discuss, in more detail, possible types of these transitions. In order to examine possible 
types of valence and insulator-metal transitions induced by pressure we have chose three 
different values of G (G = 1, 2.5, 5) that represent three typical regimes of the spin-one-half 
phase diagram. For G = 5 the valence transition can be constructed immediately. Indeed 
the ground states in this region are only the segregated configurations and thus they can be 
used directly in numerical calculations instead of the full set of /-electron configurations. 
This allows to perform numerical calculations on large systems and practically fully exclude 
finite size effects. Results of numerical calculations are shown in Fig. 4. It is seen that in 
the strong coupling limit the spin-one-half FKM exhibits a discontinuous valence transition 
from an integer valence ground state nj = 1 into another integer valence ground state 
rif = at Ef = E c ~ 1.273. The inset in Fig. 4 shows that the discontinuous valence 
transition is accompanied by a discontinuous insulator-metal transition since the energy gap 
A vanishes discontinuously at Ef = E c . Thus we can conclude that the spin-one-half FKM 
in the pressure induced case can describe the discontinuous insulator-metal transitions from 
an integer- valence state {rif = 1) into another integer- valence state (nf = 0). 

The situation for G = 2.5 is slightly complicated. Although we know that the ground 
states for 2.45 < G < 2.75 (and rif < 1/2) are configurations of a type wi,&cw e the number 
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of configurations belonging to this class is still too large for numerical calculations on 
large lattices and thus it should be further reduced. For this reason we have performed 
numerical calculations for all finite (even) clusters up to 48 sites at selected value G = 2.5. 
We have found that for any Nf > only one type of configurations and namely mixtures 
{11001100 . . . 1100}&zw e are the ground states at G = 2.5. This allows us to avoid technical 
difficulties associated with a large number of configurations and consequently to study very 
large systems (L ~ 1200). The valence transition obtained for this set of configurations (and 
Wh configurations for n/ > 1/2) is shown in Fig. 5. It is seen that the valence transition 
for G = 2.5 is very steep. A more detail analysis showed, however, that the transition 
from ra/ = lton/ = 0is not a discontinuous integer valence transition but it consists 
of several discontinuous intermediate valence transitions. Particularly, there are several 
discontinuous insulator- insulator transitions from n/ = lton/ = l/2 and a discontinuous 
insulator-metal transition from nf = 1/2 to n/ = 0. To exclude the possibility that this 
structure is a consequence of the finite size of clusters used in numerical calculations we 
have performed calculations on several large clusters (L = 400, 800, 1200) but no finite size 
effects on the valence transition were observed. 

The most complicated situation is for small and intermediate values of G. For example, 
at selected value of G = 1 the ground states are the insulating configurations Wh (for n/ > 
0.16) and the metallic configurations w a Szw e (for nf < 0.16). However, for a construction 
of the valence transition on large lattices this knowledge is too general and thus for n/ < n c 
it is again necessary to perform an additional study of the model in order to determine 
the explicit type of ground-state configurations at selected value of G. The exhaustive 
numerical studies that we have performed for n/ < n c on finite lattices up to 60 sites 
shoved that only configurations of the type w\ = {1010 rai 1010 n2 . . . 0„ 2 1010 ni 101}&w e 
or configurations of the type u>2 = {10 ni 1010 n2 1010 n3 . . . n3 1010 n2 1010 ni l}&u? e can be 
the ground states at G = 1. This fact allows us to use the set of configurations w\ 
and W2 instead of the much larger set w a &w e and consequently to construct the valence 
transition on large lattices (L ~ 1200). The results of numerical calculations obtained 
for the E/-dependence of the /-electron occupation number nf and the energy gap A 
are displayed in Fig. 6. Unlike the valence transitions obtained for G = 5 and G = 2.5 
the valence transitions in the weak coupling limit are much wider and consist of several 
consecutive discontinuous and continuous valence transitions. The valence transitions for 
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Ef < E c ~ —1.05 are insulator-insulator transitions since they realize between insulating 
ground states corresponding to the most homogeneous configurations Wh- At Ef = E c 
the spin-one-half FKM undergoes a pressure induced discontinuous valence transition from 
an intermediate valence state with n/ ~ 0.19 into another intermediate- valence state with 
rij ~ 0.14 that is accompanied by a discontinuous insulator-metal transition. Above E c the 
/-electron occupation number nf changes continuously and vanishes at Ef = Eq ~ —0.8. 

It should be noted that the weak-coupling picture of valence and metal-insulator tran- 
sitions presented in this paper strongly differs from one found in our previous work || for 
the restricted set of configurations consisting of only the most homogeneous configurations 
Wh ■ For such a set of configurations the system exhibits the metal-insulator transition only 
at Ef = Eq and this transition is continuous. This fact indicated the serious deficiency 
of the spin-one-half FKM since in some rare-earth compounds (e.g. in SmBg Q) just a 
discontinuous insulator-metal transition is observed. Fortunately, a more accurate analysis 
performed in the present paper has not confirmed this deficiency of the model. Contrary, 
our numerical results showed that the spin-one-half FKM is capable of describing not only a 
discontinuous metal-insulator transition but all basic types of valence transitions observed 
experimentally in rare-earth and transitions-metal compounds. This indicates that the 
spin-one-half FKM could yield the correct physics for describing these materials. 

In summary, the extrapolation of small-cluster exact-diagonalization calculations has 
been used to examine the ground-state phase diagram of the spin-one-half FKM in the one 
dimension. A number of remarkable results have been found, (i) The phase separation in 
the spin-one-half FKM takes place for a wide range of /-electron concentrations nf and 
d-f interactions G, including G large, (ii) In the strong coupling limit {G > 4) the model 
exhibits a pressure induced discontinuous insulator-metal transition from an integer- valence 
state (nf = 1) into another integer- valence state (nf = 0). (iii) For intermediate values 
of G (G ~ 2.5) the FKM undergoes a few discontinuous intermediate-valence transitions. 
There are several discontinuous insulator- insulator transitions from nj = ltonj = l/2 
and a discontinuous insulator- metal transition from nf = 1/2 to rif = 0. (iv) In the weak 
coupling limit (G < 2) the model undergoes a few consecutive discontinuous and continuous 
intermediate-valence transitions as well as a discontinuous metal-insulator transition. 

This work was supported by the Slovak Grant Agency VEGA under grant No. 2/4177/97. 
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Figure Captions 



Fig. 1. The ground-state phase diagram of the spin-one-half FKM obtained over the 
full set of /-electron configurations. For 1/4 < nj < 3/4 numerical calculations have been 
done on the lattice with L = 28 while for nj < 1/4 and rif > 3/4 on the lattice with L = 36. 
Four different regions of stability corresponding to mixtures w a &z,w e , w^&LWe, w c &iw e and 
WdSzwf are denoted as o, x,+ and A. 

Fig. 2. The ground-state phase diagram of the spin-one-half FKM obtained on the 
extrapolated set of /-electron configurations for L = 100. Four different regions of stability 
corresponding to mixtures w a &w e , Wb&w e , w c &w e and Wd&Wf are denoted as a,b,c and 
d. The dashed line through the exact numerical points is a guide to the eye. 

Fig. 3. /independence of the energy gap A calculated for the most homogeneous 
configurations of the type a and b (see Table 1). 

Fig. 4. Dependence of the /-electron occupation number nj on the /-level position Ef 
for G = 5 and three different values of L. Inset: The behavior of n/ and A close to the 
insulator-metal transition point. 

Fig. 5. Dependence of the /-electron occupation number n/ on the /-level position Ef 
for G = 2.5 and three different values of L. Inset: The behavior of n/ and A close to the 
insulator-metal transition point. 

Fig. 6. Dependence of the /-electron occupation number nj on the /-level position Ef 
for G = 1 and three different values of L. Inset: The behavior of nf and A close to the 
insulator-metal transition point. 
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Table Captions 

Table 1: Two basic types of the most homogeneous configurations that fill up practically 
the whole MHD for L = 24. 
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